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ABSTRACT 

We searched for X-ray shadowing toward two infrared dark clouds (IRDCs) 
using the MOS detectors on XMM-Newton to learn about the Galactic distribu- 
tion of X-ray emitting plasma. IRDCs make ideal X-ray shadowing targets of 
| keV photons due to their high column densities, relatively large angular sizes, 
and known kinematic distances. Here we focus on two clouds near 30° Galactic 
longitude at distances of 2 and 5kpc from the Sun. We derive the foreground 
and background column densities of molecular and atomic gas in the direction of 
the clouds. We find that the | keV emission must be distributed throughout the 
Galactic disk. It is therefore linked to the structure of the cooler material of the 
ISM, and to the birth of stars. 

Subject headings: ISM : clouds — Galaxy : structure — plasmas — X-rays : diffuse 
background — X-rays : ISM 



1. Introduction 

The diffuse X-ray background of the Milky Way has been studied for over 30 years and 
yet relatively little is known about the distribution of its emissi on in the Galactic plane . 



Most diffuse Galactic emission arises from ~ 1 — 3 x 10 6 K plasma ( iMcCammon et al.lll983l ). 
but the origin and distribution of this plasma are still debated. Two factors are responsible 
for the difficulty in determining the distribution of this plasma: confusion and absorption. 
Confusion arises because one cannot determine the origin of an individual X-ray. Even 
though most X-rays contributing to the diffuse background of the Milky Way are thermal in 
origin and the spectra of this plasma show strong line emission, the available X-ray detector 



institute for Astrophysical Research, Department of Astronomy, 725 Commonwealth Ave., Boston Uni- 
versity, Boston MA 02215, USA. 

2 Current address: Laboratoire d'Astrophysique de Marseille, 38 rue F. Joliot-Curie, 13388 Marseille 
Cedex 13, France 

3 NASA Goddard Space Flight Center, Code 662, Greenbelt, MD 20771, USA. 



-2- 



energy resolution is insufficient to provide velocity (and therefore distance) information. 
Thus it is very difficult to determine the true distribution of X-ray emitting plasma along a 
given line of sight. 

Shadowing experiments are the only way to determine the distribution of hot plasma 
in the plane of the Milky Way. To date, they have been used primarily in the \ keV band 
to study t he distribution of 0.1 keV plasma within the Local Hot Bubble and in the lower 

[2 keV and 



halo (e.g.. ISnowden et al.ll2000[ ). A few shadowing experiments at hig her energy 



1.5 keV) using RO SAT data (e.g.. iPark et al.lll997t lAfmy et al.l I200CH ) have shown evidence 
for the existence of an extensive distribution of hot plasma well within the solar circle. 
This emission has been linked to the Milky Way X-ray bul ge which has a scale height of 
~ 1.8 kpc and a radial extent of ~ 5 kpc (ISnowden et al.lll997l ). but there are likely additional 
components distributed throughout the Galactic disk. 

Most studies of the X-ray background have focused on the Local Hot Bubble (LHB). The 
LHB is a regi on of about 100 pc in radius, surroundi ng the Sun filled with a hot, rarefied 10 6 K 



plasma (e.g.. ISnowden et al.lll990l ). Beginning with ISnowden. McCammon. fc Verterl ( 



many groups have effectively mapped the distribution of this gas (e.g. ISmith et al 



993), 



2005; 



Galeazzi et al.l 120071 : iHenlev. Shelton. fc Kunt d 120071 : ISmith et all 120071 : iHenlev fc Shelton 
20081 ). These studies have shown the utility of X-ray shadowing in mapping a large region 
of the sky at higher Galactic latitudes. Here we utilize similar techniques but extend the 
analysis to larger distances in the Galactic plane using dense molecular clouds. 

Infrared Dark Clouds (IRDCs) are very dense molecular clouds that appear dark at 
infrared wavelengths. They are seen in absorption against the mid-IR background. IRDCs 
were first identified as a significant population in the Midcourse Space Experiment (MSX) 
mid-infrared data (jEgan et al.lll998l ). and are a ubiquitous f eature of the higher r esolution 
Galactic Legacy Infrared Survey Extraordina ire (GLIMPSE: Benjamin et al.l 120031 ) . IRDCs 
have very high column densi ties of 10 
as high as 10 23 - 10 25 cm" 2 fegan et al.lll998 



22 cm 2 (iSimon et al. 



Carey et al 



2006 . hereafter S06), and possibly 
19981 ) . They exist at a wide range 



of distances and are most easily seen in the plane of the Galaxy (S06). 

S06 found distances to 313 IRDCs using the Galactic Ring Survey of 13 CO emission 
(GRS; I Jackson et al.l 120061 ) to establish a morphological match between th e 13 CO gas and 
the mid- infrared extinction. Using the rotation curve of I Clemens! (119851 ). S06 converted 
the velocity of the associated 13 CO gas into a distance to the IRDC. In the inner Galaxy 
(the focus of the present study), every velocity has two possible distance solutions. IRDCs, 
however, are seen in absorption against the Galactic plane and therefore the near distance 
can be assumed. 
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IRDCs make ideal X-ray shadowing candidates. Their high column densities make 
them effective absorbers of background X-ray photons; X-ray emission from the direction 
of the IRDC must originate in the foreground. Thus, using IRDC absorption, one can 
separate foreground and background X-ray components. The distribution of IRDCs also 
spans a wide range of distances (see S06. 1 Jackson et al.ll2008l ). With enough X-ray absorption 
measurements of IRDCs along different lines of sight, an accurate map of the hot ISM may 
be created. 



2. IRDC Sample 

The present study includes two clouds from the catalog of S06 that have high column 
densities, large angular sizes, and are closely grouped on the sky. Our goal is to select clouds 
that would cau se complete absorptio n of X-ray photons originating behind the cloud. The 



models given in lSnowden et al.l ( 11994] ) predict that a total column density of 1 x 10 22 cm -2 will 
absorb ~ 98% of incident photons at | keV. Column densities of 10 22 are frequently found 
for IRDCs (see S06). To ensure sufficient counts in both the on- and off-cloud directions, 
the target IRDCs should cover approximately half the XMM-Newton detector. Our ideal 
situation is one in which there are multiple large clouds closely grouped on the sky but at 
various distances. The close grouping allows one to disentangle the emission components 
along a given line of sight. 

The properties of the two IRDCs observed by XMM-Newton are summarized in Tabled] 
Listed are the IRDC name, the Galactic longitude and latitude, the angular size, and the 
molecular column density. All parameters in Table [JJ are reproduced from S06. The two 
target clouds have high column densities, large angular sizes, and are near / = 30°. 

The mid- infrared emission of the two IRDC fields is shown in Figure [U Th e background 



image in this figure is from the GLIMPSE 8//m survey ( jBenjamin et al.ll2003l ) and the con- 
tours are 13 CO integrated intensity from the GRS. We have created the integrated intensity 
map using the LSR velocity and line width from the S06 catalog. Tick-marks on the contours 
indicate the direction of decreasing 13 CO integrated intensity. The large circle in this figure 
is approximately the XMM-Newton field of view (~ 30'). Immediately evident in this figure 
is the excellent morphological match between the GRS 13 CO contours and the strong mid- 
infrared extinction. The molecular clouds associated with the IRDCs take up a significant 
fraction of the XMM-Newton field of view. IRDC G28.37+0.07 is located at the center of a 
large, dense molecular cloud which is roughly ellipsoidal in projection. The two clouds are 
at very different distances, G28.37+0.07 lies ~ 5 kpc from the Sun and G36.67— 0.11 lies 
~ 2 kpc from the Sun (S06). 
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Observations and Data Reduction 



Our clouds were observed by XMM-Newton on 2006 April 11 for 25 ks each (IRDC 
G28.37+0.07 with XMM-Newton ObsID 0302970301 and IRDC G36.67-0.ll with ObsID 
0302970201). 

One of us (SLS) has developed a powerful suite of PERL and FORTRAN scr ipts to 
analyze XMM-Newton extended source data called XMM- ESAS (ISnowden et al.ll2008 ) . This 
soft ware is based on the backg round modeling described in lSnowden. Collier fc Kuntzl (120041 ) 
and iKuntz fc Snowdenl (120081 ). The software currently only operates on the MOS detector 
data and subsequent to this analysis has been incorporated into the XMM-Newton mission 
Standard Analysis Software (SAS) package. 

To remove intervals of high soft proton contamination, we must first temporally filter 
and clean the events file. The signature of soft proton contamination is a fluctuating light 
curve. The XMM-ESAS software creates a count rate histogram and fits a Gaussian to 
the peak of the histogram. In observations with minimal contamination, this histogram 
would appear Gaussian; a non-Gaussian shape to the distribution is also indicative of soft- 
proton contamination. The software defines count-rate intervals within 1.5a of the median 
count rate as acceptable, and filters the events file to remove intervals that do not satisfy 
this criterion. While effective at removing intermittent contamination, this process does 
not remove consistent low-level contamination that is also likely present. The light-curve 
for G28.37+0.07 is quite steady and almost all the data are useable. The light-curve of 
G36.67— 0.11, however, indicates a high level of soft proton contamination and only ~ 20% 
of the observation time is usable. 



4. Estimation of Column Density 

An estimate of the total column density to each IRDC can drastically improve the 
reliability and accuracy of the X-ray spectral fitting. Above 0.5 keV, the most significant 
source of X-ray absorption is metals. With normal abundances, the total X-ray opacity 
should be proportional to iV(H) + 2A^(H 2 ). Thus, we estimate the molecular column density 
using the 13 CO emission from the GRS, and es timate the atom ic column density using the 



21 cm Hi VLA Galactic Plane Survey (VGPS: IStil et all 120061 1 . 



Our spectral fitting procedure treats the on- and off-cloud spectra separately (as do all 
X-ray shadowing experiments). This is the only way one can separate components of the 
X-ray background that originate in front of the cloud from those that originate behind it. 
One must therefore compute on- and off-cloud column densities separately. We do this by 
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first defining the on- and off-cloud regions. 

One of us (LDA) has developed an astronomical software program in IDL§ This soft- 
ware allows one to define regions of interest in an image using a threshold selection tool. 
The user can vary the threshold upwards or downwards to define smaller or larger regions. 
The output of this thresholding selection, the x and y positions of the data points on the 
threshold boundary in sky coordinates, can then be converted to detector coordinates within 
the software, as is appropriate for analysis in the XMM-Newton Science Analysis Software 
(SA£@). This conversion employs the SAS program edet2sky repeatedly for each coordinate 
pair to convert from a region in sky coordinates to a region in detector coordinates. 

Using our IDL software, we define the on-cloud region as the sky positions possessing 
both high mid-infrared extinction and high 13 CO column density. Our goal is to create a re- 
gion large enough to have sufficient counts in the on-cloud region for good statistics. Ideally, 
the on- and off-cloud regions would each cover half the XMM-Newton field of view. The 
off-cloud region encompasses all locations in the XMM-Newton field of view not including 
the on-cloud region. The on-cloud region roughly corresponds to the largest contour in 
Figure [TJ 



4.1. Molecular Column Density 

We must separate the total molecular column density into a component foreground to 
the IRDC and a component beyond the IRDC We therefore must estimate the distances to 
the molecular gas clumps that contribute to the column density in the XMM-Newton field 
of view. In the first Galactic quadrant (where the target IRDCs lie), each radial velocity 
has two possible distance solutions, a "near" and a "far" distance. This problem is known 
as the kinematic distance ambiguity (KDA). Gas along the same line of sight as an IRDC 
could lie either in front of or behind the cloud. This problem is illustrated in Figure [2] for 
IRDC G28.37+0.07. Figure [2] shows the loci of distances fr om the Sun cor responding to a 



particular velocity, assuming the Galactic rotation curve of IClemensI (119851 ). The two- fold 
distance degeneracy is clearly shown. Because IRDCs are seen in absorption, we assume 
they lie at the near distance. The cloud velocity of 78.6 kms -1 places G28.37+0.07 at a 
distance of 4.8 kpc from the Sun. Any gas with a velocity greater than 78.6 kms -1 must 
lie beyond the cloud. Molecular gas at velocities less than 78.6 kms -1 , however, could lie 



1 Download at http://www.bu.edu/iar/kang/ 
2 http://xmm.esac.esa.int/sas/ 
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at either the near or the far distance. 



Hi self-absorption (Hi SA) can remove the distance ambiguity of the molecular gas 
clumps in the XMM-Newton field of view. This technique relies on cold foreground H I ab- 
sorbing the emission of warm background H I at th e same velocity. Warm H I is ubiquitous 
in our Galaxy and emits at all allowed velocities. iLitzt. Burton fc Banial (119811 ) hypothe- 
sized that molecular clouds mu st contain H I, a result which has been confirmed in nu merous 
subsequent observations (e.g.. iKuchar fc Banial 1 19931 ; IWilliams fc Maddalenal Il99q ). This 
small population of cold neutral Hi atoms inside molecular clouds is maintained by inter- 
actions with cosmic rays. The H I inside molecular clouds at the near distance will produce 
an absorption signal because there is ample warm Hi emitting at the same velocity at the 
far distance. Any cloud at the far distance should not show Hi SA because there is no 
background H I emitting at the same velocity. H I SA is therefore a general technique for 



finding distances to molecular clouds. This technique was shown b y I Jackson et a 



be effective at resolving the KDA for a dense molecul ar cloud, by iBusfield et al 



resolving the KDA for young stellar objects, and by lAnderson fc Banial (120091) for resolv 



(2002|) to 



(J2006|) for 



i ng th e KDA for the molecular gas associated with H II regions, and by iRoman-Duval et al. 
(120091 ) for resolving the KDA for GRS molecular clouds. 



We utilize the procedure outlined in ISimon et al.l (120011 ) to convert from GRS 13 CO 
emission line parameters into an estimate of column density. This procedure assumes a 
standard CO excitation temperature of 10 K. The conversion equation is: 



iV(H 2 ) = 5.2 x 10 20 TAV [cm 



(1) 



where T is the main beam line intensity in degrees Kelvin and AV is the line width in 
kms -1 . The column density found using this procedure is accurate to within a factor of ~ 2 
due to uncertainties in the conversion from 13 CO to H2. 

To find the molecular column densities for both on- and off-cloud directions, we calcu- 
late the column density of all significant 13 CO molecular clumps in the XMM-Newton field of 
view by first fitting Gaussians to the average on- and off- cloud GRS spectra. We transform 
the line parameters derived from these Gaussians into column densities using Equation [TJ 
The application of Equation [1] computes the total column density along each line of sight, but 
gives no information about whether the components lie at the near or the far distance. By 
examining the average 13 CO and Hi spectrum of each molecular clump in the XMM-Newton 
field of view for H I SA, we determine whether the near or the far distance is appropriate for 
each molecular clump. We assign the column densities of the molecular clumps found to lie 
in front of the target cloud to the foreground component, and the column densities of the 
molecular clumps found to lie behind to the target cloud to the background component. 
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4.2. Atomic Column Density 

We use a similar procedure to calculate the Hi column density using the VGPS. The 
conversion to Hi column density employs the equation: 



iV(HI) = 1.82 x 10 18 T s J t(v) dv [cm" 



(2) 



where T s is the Hi spin temperature in degrees Kelvin, and v is given in kms -1 . For Hi, 
the spin temperature is much less well constrained compared to the excitation temperature 
of 13 CO. A good averag e value between the hot and cool Hi components is 150 K (see 



Dickey fc Lockmanlll990l . and references therein). Equation requires an estimate of the 



optical depth for the H I line using the equation: 

where thi is the optical depth of H I, Tl is the H I line intensity, T s is the H I spin temperature, 
and Tbg is the intensity from background sources. The largest source of background emission 
is the cosmic microwave background; because it produces a negligible effect at 21 cm, we 
ignore its contribution here. Assuming a spin temperature of 40 K lowers our calculated 
column densities by a factor of ~ 2. 

For H I there is less information about whether the gas originates in front of or behind the 
cloud compared to CO. As we did for 13 CO, we compute average on- and off-cloud spectra. A 
typical H I spectrum does not show distinct, clean lines. We therefore compute the area under 
the curve directly, instead of fitting Gaussians, and convert this to H I column density using 
Equations [2] and El We assume that the H I is uniformly distributed along the line of sight to 
a distance of 15 kpc from the Sun and split the column density distribution into foreground 
and background components, weighted by the distance to the cloud. For example, we assume 
G28.37+0.07 (at a distance of 5 kpc) has 33% of the total Hi along the line of sight in the 
foreground component, and 67% in the background component. Our values for the total Hi 
column de nsity are in very good agreement with thos e in the the Leiden/Argent ine/Bonn 
Hi survey ( Kalberla et al. 2005^ and the Hi survey of Dickey Sz Lockman ( 1990 ) as found 



usm, 



g the HEASARC online tool@ 



Figure [3] shows the average on- (solid curves) and off-cloud (dashed curves) spectra of 
13 CO (black curves) and Hi (gray curves) for Figure El For G28.37+0.07, Hi absorption is 
seen at the velocity of all 13 CO components with velocities less than the cloud velocity of 



3 http://heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 
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78.6 km s -1 (this is clearer when analyzing each clump individually). Therefore, we assign 
all 13 C0 emission not associated with G28.37+0.07 that has velocities less than 76.6 kms -1 
to the near component. We assign all 13 C0 features with velocities greater than 78.6 kms -1 
to the far component. Also evident in Figure [3] is the strong Hi absorption at the cloud 
velocity, which shows clearly that the cloud is at the near distance. The H I emission is less 
strongly peaked than the 13 CO emission. 

The results of our column density modeling are summarized in Table [2J This table 
gives the column density of atomic hydrogen, iV(H), molecular hydrogen 7V(H 2 ), and the 
total column density, A Tot ai = A(H) + 2A(H 2 ) in units of 10 22 cm~ 2 for G28. 37+0.07 and 
G36.67— 0.11. The contributions to the column density are given for the foreground and 
background components for the on- and off-cloud directions. The background component 
includes the contributions from the IRDCs themselves. 



5. Spectral Analysis 

Using our cleaned events files, we extract the X-ray spectra for both the on- and off- 
cloud regions in preparation for spectral fitting. Our on- and off-cloud regions are defined as 
before and correspond to high infrared extinction and high 13 CO column density. We extract 
spectra from both MOS detectors and create auxiliary response files (ARFs), redistribution 
matrix files (RMFs), and background spectra. We increase the signal to noise by binning 
these spectra so the re are at least 20 counts per spectral channel. We use XSPEOj Vll.3.2 



( lArnaud et al.l 119961 ) for the spectral analysis. 

For the spectral fits, we assume the emission can be described as: 

I x = SP + I LHB + I F + I B + (I R + AGN) , (4) 

where Ix is the measured X-ray flux, SP is the contribution from soft proton contamination, 
Ilhb is the emission from the LHB, Ip is the Galactic X-ray emission foreground to the 
cloud, Ib is the Galactic X-ray emission background to the cloud, Ir is the emission from 
the Galactic ridge, and AGN is the X-ray background from unresolved active galactic nuclei. 
There are three column density parameters in our models, one associated with Ip, one with 
Ib, and one with both Ir and AGN. The fitted column density associated with Ib is poorly 
constrained (see below), which is why we do not treat the column density associated with 
Ir and AGN as the sum of that affecting Ip and Ib- Below we describe in more detail the 



http : / /xspec . gsfc . nasa.gov/docs /xanadu /xspec 
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various model components we use to model the X-ray background and to separate foreground 
from background emission. 

There are two main sources of contamination in our spectra. First, there is residual 
soft proton contamination not removed through time filtering (SP in Equation H]). This 
component should have a power law distribution of energies, so we model its emission using 
the XSPEC model pow/t£. This XSPEC component models a power-law noise distribution 
not folded through the instrument response. Secondly, there are two fluorescent instrumental 
lines of Al Ka and Si Ka near 1.6 keV. We assume a Gaussian shape for these lines and fit 
them explicitly. 



Our models include three components for the X-ray background, as in ISnowden et al. 



(2000). Along each line of sight, first there is ~ 0.1 keV emission from the LHB (ISmith fc Cox 
200lh . I lhb- We fit the the LHB emission with the XSPEC model apec JSmith et allboOlh . 
which represents collisionally-ionized diffuse plasma. Next, we assume that there are two 
components of the X-Ray background emission in the Galactic plane: one soft (< 1 keV) 
and one hard (few keV), Jr. For these components we use the XSPEC model wabs x apec, 
which represents collisionally-ionized diffuse plasma attenuated by material along the line 
of sight. The soft component itself can be divided into two sources: emission in front of 
the IRDC, If, and emission from beyond the IRDC, Ib- Finally, there is extragalactic 
emission from unresolved AGN in the field, AGN. We model the contribution from AGN 
using a p ower law model w ith absorption, wabs x pow, where the power law index is set to 
a = 1.46 (IChen et al .11199 71 ) . This is expected to be the dominant component above 1 keV 



iLumb et al 



2002 



We show the model components and the column densities affecting them graphically in 
Figure H] for G28.37+0.07. In this figure, molecular gas is displayed in black and atomic gas 
is displayed in gray. The white wedge is an exaggerated representation of the XMM-Newton 
field of view. Also shown on the figure are the location of IRDC G28.37+0.07, the LHB, 
and the direction to the Galactic center (GC). The LHB size in the figure is exaggerated 
for clarity. Figure H] shows that the AGN X-ray emission is affected by the total Galactic 
column density. 

During the spectral modeling, we fit the spectra for each cloud simultaneously to improve 
the fit quality. We fix the abundance to Solar and use the Solar relative metal abundance of 
Anders fc Grevessd (119891 ). We have five spectra for each cloud: on- and off-cloud spectra 
for both MOS detectors, as well as a ROSAT All Sky Survey (RASS; ISnowden et al.lll997l ) 



1 The notation for this model is that of XSPEC Vll; in more recent versions the "/b" has been dropped. 



See http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/backgroundmodel. htaiil 



- 10 - 



spectrum derived from the HEASARC X-ray Background Toolo The RASS spectrum helps 
constrain the cosmic background at lower energies where the XMM-Newton detectors are 
less sensitive. We extract the RASS spectra from regions centered on our clouds of radius 
0.5° and normalize the spectra to 1 arcminute 2 . We use the same RASS spectrum for the 
on- and off-cloud regions, but scale the emission by the relative area of these regions. We 
link the normalizations and temperatures of each spectral component for all five spectra to 
constrain the fits. We also link the temperature of the soft foreground component to that of 
the soft background component. 

The results of our spectral fitting procedure are shown in Figure|5]for the MOS1 detector. 
In Figure the top curve is the on-cloud spectrum, and the lower curve is the off-cloud 
spectrum. The bottom curve is the RASS spectrum. The two prominent emission lines are 
the fluorescent instrumental lines of Al Ka and Si Ka. The left panel of Figure shows the 
spectra and model for G28.37+0.07, and the right panel shows the same for G36.67— 0.11. 
It is evident from this figure that our model adequately fits the data. Also apparent from 
Figure [5] are the decreased counts for G36.67— 0.11 due to the soft proton contamination. 
This results in greater uncertainty in the spectral fits. 

For G28.37+0.07, we run our fits in two trials, once with the column densities fixed to 
their calculated values (see §1]) and once with the column densities allowed to vary (hereafter 
the "fixed" and "free" trials). We allow the column densities to vary up to 100% from 
our calculated values to account for uncertainties in converting line parameters to column 
densities. Due to insufficient counts, we are unable to allow the column densities to vary for 
G36.67— 0.11, and thus we run the spectral fits only with fixed values. 

5.1. Spectral Fit Results 

The results of our spectral fits are shown in Figure [6] and listed in Table [31 Figure [6] 
shows the best fit model components for G28.37+0.07 (left panels) and G36.67— 0.11 (right 
panels). The top panels show the models for the on-cloud directions, and the bottom panels 
show the same for the off-cloud directions. The models for G28.37+0.07 are from the free 
trial. Table [3] lists the fitted temperature and emission measure (EM = J n 2 e dl) for our four 
spectral model components, the reduced x 2 value, and the number of degrees of freedom in 
the fit. Errors in the table represent 90% confidence levels. Below we discuss the model 
components individually. 



5 http:/ /rosat. gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg. pi 
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5.1.1. Column Densities 

When allowed to vary, the column densities for G28.37+0.07 fit to similar, but generally 
lower, values compared to our calculated values (see For the on-cloud spectra, the 

column density fits to 1.21qq2 x 10 22 cm" 2 for the foreground component; we estimated 
1.3 x 10 22 cm -2 . The column density background to the cloud is unconstrained when allowed 
to vary as its normalization fits to zero (see §5.1.4l below). The total column density affecting 
the Galactic ridge and AGN emission fits to 3.9i£j x 10 22 cm for the on-cloud spectra; 
we estimated 5.2 x 10 22 cm" 2 . Thus, since the total column density affecting the Galactic 
Ridge and AGN component is the sum of that foreground and background to the cloud, we 
find 2.6 x 10 22 for the background component. For the off-cloud spectra, the column density 
fits to l.liaol x 1q22 cm ~ 2 foreground to the cloud and to 2.9^2 x 10 22 cm" 2 for the total 
column density; we estimated 0.9 x 10 22 cm -2 and 4.2 x 10 22 cm -2 , respectively. 

Perhaps the most important column density parameter - the difference in column den- 
sity between the on- and off-cloud regions - is the same in the fixed and free trials for 
G28.37+0.07: 1.0 x 10 22 cm -2 , although with a formal error derived from the fit of about 
50%. The other discrepancies are discouraging, although the differences have little affect on 
the other model parameters. The largest effect can be seen in the emission measure of the 
foreground soft component (see below). Despite our care taken in estimating the column 
densities, errors in the conversion from 13 CO line parameters to hydrogen column density 
are large and currently unavoidable. Longer or more sensitive observations are required to 
more accurately disentangle the competing effects of absorption from the Galactic column 
density and emission from the X-ray background. 



5.1.2. LHB Component 



The fixed and free column density trials for G28.37+0.07 both result in a LHB compo- 
nent temperature kT = 0.11 keV (logT/K = 6.11). For G36.67— 0.11, the LHB component 
fits to a temperature of 0.089 keV(logT/K = 6.02. These temperatures for t he LHB compo- 



nent are in good agreement with that found by many previous stu dies (see ISnowden et al. 



19981 . 120001 ; iKuntz Snowdenll2000i ; iHenley. Shelton. fc Kuntzll2007l ) . The emission measure 
for the LHB component is less than what has been found by previous authors. We calculate 
an emission measure of 5.5 x 10~ 4 cm" 6 pc for the fixed trial and 5.4 x 10~ 4 cm -6 pc for 
the free trial for G28.37+0.07. For G36.67-0.ll, it is slightly higher: 6.3 x IP" 4 cm 



pc. 
10 



These values are three times le ss than that f ound by IHenley. Shelton. fc Kunta (120071 ). 
times less than that found by ISmith et al.l (120071 ). This emission measure is in large part 
determined from the RASS 1/4 keV band data, and as the directions lie in the Galactic plane 



the surface brightness, and therefore the emission measure of the RASS data are relatively 
low. 



5.1.3. Soft Foreground Component 

The soft component fits to a temperature of 0.29 keV (logT/K = 6.53) in the fixed trial 
and to 0.30 keV (logT/K = 6.54) in the free trial for G28.37+0.07. For G36.67-0.ll, the 
soft component is not very well constrained. In fact, the only c onstraint we can p l ace on 
this component is that it is less than 0.35 keV(logT/K = 6.61). Kuntz fc Snowdenl (120081 ) 
find a similar temperature for this component: 0.24 keV. 



5.1.4- Soft Background Component 

For the fixed trial of G28.37+0.07, the soft background component has little effect on 
the model. In fact, its removal has no impact on the quality of the fits, nor any impact 
on the value of the other parameters. The large column density in front of the background 
component causes significant absorption. The negligible effect this component has on our 
total model can be seen in Figure [6j For the free trial, the normalization for this compo- 
nent fits to zero. For G36.67— 0.11, the soft foreground component and soft background 
component have similar intensities. The difference between the fits seen in Figure |6] for the 
two clouds is due partly to the decreased column density to G36.67— 0.11. Perhaps a larger 
factor, however, could be the greater uncertainty in the spectral fits for G36.67— 0.11 due to 
soft proton contamination. 



5.1.5. Galactic Ridge Component 

The temperature for the warm Galactic ridge component fits to 1.61 keV (logT/K 
7.27) in the free trial and to 2.34 keV (logT/K = 7.43) in the fixed trial for G28.37+0.07. This 



i s a si gnificantly higher than what has been found by previous authors. Kuntz fc Snowden 



( 120081 ) find a temperature of 0.71 keV for this component, with rather large error bars. 
Because of the decreased counts in the observation of G36.67— 0.11, we were forced to fix 
the Galactic ridge component to the temperature found for G28.37+0.07. When allowed to 
vary, this component would fit to an extremely high temperature (> 10 keV), which we deem 
unphysical, especially considering the results for G28.37+0.07. 



Figure |6] shows that for both clouds the Galactic ridge emission is the strongest compo- 
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nent of the model between 2 and 5 keV. To determine if this component is required by the 
spectral fits, we remove the hard component from the models for G28.37+0.07. Its exclusion 
results in significantly worse fits: the reduced x 2 value increases from 1.11 to 1.41. An F-test 
reveals that the inclusion of the Galactic ridge component is significant at the 15a level. We 
conclude that this is a necessary component of our spectral fits although its temperature is 
not very well constrained. 



5.2. O VII and O VIII Emission 



Oxygen is the most abundant line-emitting element at the temperatures of the LHB. In 
a million degree plasma, oxygen is primarily in the +6 state (O vil), and the strongest lines 
detectable by XMM-Newton are near 0.5 keV. The emission from other charge states (most 
important for our purposes is O Vlll) can be used as a diagnostic to constrain the temperature 
of the LHB and verify the temperature found previously through spectral fitting. There are 
three lines of O vil that are all near 0.57 keV, and there is a Lya transition of O Vlll at 0.65 
keV. 

We measure the intensities of the O VII and O Vlll emission from the LHB by replacing 
the apec model component with a vapec model. The vapec model is identical to the apec 
model except that the atomic abundances may be individually modified. We set the oxygen 
abundance to zero and add two Gaussians representing the O VII and O Vlll emission lines 
at 0.57 keV and 0.65 keV. For G28.37+0.07, when we fix all the model parameters to their 
previously fit values (except for the vapec temperature and normalization), we find that 
the Ovil line has an intensity of 3.4l5'o photons cm~ 2 s _1 sr _1 (hereafter line units, LU). 
The vapec temperature fits to a value of 0.098 keV. When we fix the vapec temperature to 
the LHB temperature found previously, 0.11 keV, we find an Ovil intensity of 3.3^8 LU. 
We do not detect the O Vlll line in either of these trials. The same procedure applied to 
G36.67— 0.11 yields similar results. With the vapec temperature allowed to vary, the Ovil 
intensity is 1.9lJ'g LU. The vapec temperature fits to 0.085 keV. When fixed to the 0.089 
keV value found previously, the Ovil intensity is 2-l^J LU. We do not detect the Ovill 
line. 

Our values for the O vil and O Vlll lines are in rough agreement with what has been 



found by other authors. 



n_a shadowing experiment to a nearby (~ 200 pc) filament, 



Henlev. Shelton. fc Kuntzl (120071) found 3.4±g;| LU for the Ovil line and did not detect 



the Ovm line. iKuntz fc Snowdenl (120081 ) found 1.75 ±0.7 LU for the Ovil emission at (l,b) 
= (111P14, If 11). In a shadowing e xperiment of MBM 12, a high latitude cloud ~ 50 to 
~ 300 pc distant, ISmith et al.l fepQsh found 1.79 ± 0.55 LU for O VII and 2.34 ± 0.36 LU for 
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O vm. They note that this unusually strong O Vlll emission can be explained by solar wind 
charge exchange (SWCX) contamination in their spectra. The SWCX emission is caused 
by electron transitions between neutral atoms and highly ionized species in the solar wind. 
It has a potentially large effect in X-ray data, especially following a coronal mass ejection 



event. ISnowden. Collier fc Kuntzl (120041 ) showed the large effect on th e oxygen lines that 
SWCX emission can have. In a followup study of MBM 12 with Suzaku, ISmith et all (120071 ) 
found 2.34 ± 0.33 LU for O VII and 0.77 ± 0.16 L U for O vm line em i ssion Using the X-ray 
Quantum Calorimeter aboard a sounding rocket, iMcCammon et all (120021 ) found 4.8 ± 0.8 
LU for the O VII emission and 1.6 ± 0.4 LU for the O vm emission. 

While consistent with previous values, we note that there is likely some resi dual contam- 
i natio n from SWCX in our measured oxygen line intensities. Using the models in lSmith et al. 
( 120051 ). we estimate that for the EM and T found here for the LHB, we should expect O VII 
and O vm intensities of ~ 2 and < 1 LU, respectively. This indicates that the line intensities 
measured here cannot be solely due to a plasma in thermal equilibrium. 



A consensus appears to be emerging that most of the O VII emission arises from the 
LHB (with a contribution from SWCX). Strong O Vlll emission, however, cannot arise from 
the LHB as this would imply a temperature inconsistent with o ther determinations. This 



emission must arise from SWCX (ISnowden. Collier Kuntz 120041 ) or from a cool component 



distributed throughout the Galaxy ( IKuntz fc Snowdenll2008l ). likely with contributions from 
both sources. Our results are consistent with this interpretation. That we do not detect the 
O vm line is not surprising given the large absorbing column density and low expected line 
strength. 



6. Image Analysis 

We produce exposure corrected, background subtracted X-ray images to verify the re- 
sults from the spectral analysis, and to search for a morphological match between the molec- 
ular emission and an X-ray decrement. Using the XMM-ESAS software, we create these 
images in three energy bands: 0.35-1.25 keV, 1.25-2.00 keV, and 2-10 keV. We model two 
sources of background in our images: the particle background and residual soft proton con- 
tamination. These background components must be removed as cleanly and completely as 
possible to ensure accurate results. 



Using t he filter-wheel-closed da ta, the XMM-ESAS software models the particle back- 
ground (see lKuntz &: Snowdenll2008l ). These data are dominated by the instrumental back- 
ground because the chips were not exposed to the sky. We do not include blank sky data in 
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the modeling of either the particle background nor the soft proton contamination. The blank 
sky files may suffer from residual soft proton contamination and solar wind charge exchange. 
Additionally, the blank sky files contain the cosmic background which we are attempting to 
observe. Using these data would not only add error to our models, but would also remove 
the signal of interest. 

The XMM-ESAS softwa re also allows for th e creat ion of a model image of any residual 



soft proton contamination. iKuntz fc Snowdenl (120081 ) characterize the reasonable ranges 



of the variations in both the spectrum and in the spatial distribution over the detectors. 
XMM-ESAS incorporates this characterization in the model soft proton image. The energy 
distribution of the soft proton contamination is determined from the power law index and 
normalization found in the spectral fits. For this input, we use the values found in the 
spectral fitting. 

We subtract the model background and soft proton images, combine the data from the 
two MOS detectors and smooth the resultant image. The result of these operations is an 
exposure corrected, background subtracted image. The strongest X-ray shadow should be 
present in the 0.35-1.25 keV image. We show in Figure[7]the exposure corrected, background 
subtracted image for the 0.35-1.25 keV band of G28.37+0.07 (left panel) G36.67-0.ll (right 
panel). In the images of both clouds there is some indication of a shadow, but it not the 
clearly defined shadow one would expect, nor does it have the same depth one would expect. 
Furthermore, the same shadow that appears in the 0.35-1.25 band also appears in the 2.0-8.0 
keV band. As there should not be strong absorption in the 2.0-8.0 keV band, this indicates 
that the "shadow" is a product of the detector, and not of the cloud itself. We conclude that 
the imaging analysis supports the spectral analysis and finds no X-ray shadow of the cooler 
component to the X-ray background. 



7. Summary 



We have conducted X-ray shadowing experiments on two infrared dark clouds near 
/ = 30° at distances of 2 and 5 kpc from the Sun to determine the spatial distribution of the 
X-ray background. We used Hi and 13 CO data to calculate the total atomic and molecular 
column densities for both the on- and off-cloud directions. Both clouds have very high 
proton column densities of > 4 x 10 22 cm~ 2 . This high column density should absorb nearly 
all the soft background X-ray flux. We find that the cool, diffuse X-ray background must 
originate foreground to the clouds, within a few kpc of the Sun. Our results are further 
evidence that X-ray emitting plasma is d istributed throughout the disk of our Galaxy (see 
Park et al.lll997l : IKuntz fc Snowdenl l2008h . 
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We found that the X-ray background is best fit with a three component model with 
contributions from the Local Hot Bubble (LHB), a soft component, and a hot component 
from the Galactic ridge. The LHB emission is best fit with a value of ~ 0.1 keV (log T/K = 
6.06) for both clouds, in agreement with what has been found in previous studies. The soft 
component is best fit with a temperature of ~ 0.3 keV (log T/K = 6.54). This is the dominant 
source of emission between 0.7 and 1.0 keV. The Galactic ridge component is best fit with a 
temperature of 2 keV. This is the dominant source of emission in our spectral models between 
2 and 5 keV. Our spectral fits show that most significant emission below ~ 0.7 keV can be 
attributed to the LHB (and/or to SWCX). Image analysis including accurate modeling of 
any background components failed to reveal an X-ray shadow, supporting the results of the 
spectral modeling. 

We detect the Ovn line for both clouds at levels of ~ 3 photons cm" 2 s" 1 sr" 1 . This 
intensity is ~ 50% higher than what one would expect from a thermal plasma with the 
emission measure and temperature of the LHB, indicating some level of contamination from 
solar wind charge exchange. These values are roughly in agreement with what has been 
found in previous studies. We do not detect the O Vlll line for either cloud. 

This research was based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by ESA Member States 
and NASA, and was supported by NASA XMM-Newton Guest Observer grants including 
NNX06AG73G and NNG05GP68G. We make use of molecular line data from the Boston 
University-FCRAO Galactic Ring Survey (GRS). The GRS is a joint project of Boston 
University and the Five College Radio Astronomy Observatory (FCRAO), funded by the 
National Science Foundation under grants AST-9800334, AST-0098562, AST-0100793, AST- 
0228993, & AST-0507657. We also make use of data from the VLA Galactic plane survey 
(VGPS). The National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. We would 
like to thank the referee for their careful reading of this paper and their helpful comments 
which greatly improved its clarity. 
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Fig. 1. — IRDCs shown in absorption against GLIMPSE 8/iin emission. The left panel shows 
G28.37+0.07 and the right panel shows G36.67-0.ll. The contours are 13 CO integrated 
intensity from the GRS. Tick-marks on the contours point downhill, toward decreasing 13 CO 
integrated intensity. The large black circles represent the XMM-Newton field of view (~ 30'). 
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Fig. 2. — Velocity versus distance plot for G28.37+0.07 showing the kinematic distance 
ambiguity. This cloud lies at the near distance of 4.82 kpc. 
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Fig. 4. — Components of the X-ray background attenuated by the Galactic column density 
for G28.37+0.07. The column density Molecular gas from the GRS is shown black while 
atomic gas from the VGPS is shown in gray. The Local Hot Bubble (labeled LHB), the 
direction of our observation (white wedge) and the direction towards the Galactic center are 
also shown. The size of the wedge and the LHB are slightly exaggerated for clarity. IRDC 
G28.37+0.07 is associated with a large dense molecular cloud 5 kpc from the Sun.) 
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Fig. 5. — Spectral data and models for G28.37+0.07 (left panel, column density free fit) 
and G36.67— 0.11 (right panel) from the MOS1 detector. Shown are the on-cloud (black), 
off-cloud (dark gray), and RASS (light gray) spectra. The bottom section of each panel 
shows the fit residuals. 
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Fig. 6. — Model spectral components from the XSPEC fit to the on-cloud (top panels) 
and off-cloud (bottom panels) spectra for G28.37+0.07 (left panels) and G36.67— 0.11 (right 
panels) from the MOS1 detector. For G28.37+0.07, the models are from the trial with 
column density free. 
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Fig. 7. — Smoothed XMM-Newton images in the energy range 0.35 keV to 1.25 keV of 
G28.37+0.07 (left panel) and G36.67— 0.11 (right panel). The contours are the same GRS 
contours shown in Figure [TJ There is no obvious X-ray decrement at the location of the 
molecular clouds. 
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Table I. Infrared Dark Cloud Properties 



Name 


1 


b 


D 


Size 


N(H 2 ) 




deg. 


deg. 


kpc 


sq. arcmin 


10 22 cm" 2 


G28.37+0.07 


28.37 


+0.07 


5.0 


42 


21 


G36.67-0.ll 


36.67 


-0.11 


2.0 


19 


19 



Table 2. Cloud Column Densities''" 



Name 


AT(H) 


N(H 2 ) 


AT(Total) 




FG 


BG 


FG 


BG 


FG 


BG 


G28.37+0.07 














On-cloud 


0.5 


1.0 


0.4 


2.1 


1.3 


5.2 


Off-cloud 


0.5 


1.0 


0.2 


1.6 


0.9 


4.2 


G36.67-0.ll 














On-cloud 


0.5 


1.0 


0.1 


1.4 


0.7 


3.8 


Off-cloud 


0.5 


1.0 


0.1 


0.9 


0.7 


2.8 


^ All column 


densities 


arc in units 


ol 10 


22 cm~ 2 . 


"FG" and 


"BG" 



represent foreground and background to the cloud. 
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Table 3. Spectral Model Parameters 



Parameter 



G28.37+0.07 



Column Fixed 



Value 



Error 



Column Free 



Value 



Error 



G36.67-0.ll 



Column Fixed 



Value 



Error 



Column Density X 10 22 cm 
On-cloud, Foreground 
On-cloud, Background 
Off-cloud, Foreground 
Off-cloud, Background 

LHB 

kT(keV) 

EM(cm~ 6 pc) 

Cool Foreground 
kT(keV) 

EM(cm~ 6 pc) 
Cool Background 
kT(keV) 

EM(cm~ 6 pc) 
Galactic Ridge 
kT(keV) 

EM(cm~ 6 pc) 

X 2 
dof 



1.3 
5.2 
0.9 
4.2 

0.11 
5.5E-4 

0.29 
4.6E-2 

0.29 
1.9E-3 

1.61 

2.2E-2 
1373.5 
1236 



+3.7E-3 
-9.3E-3 
+ 7.6E-5 
-7.0E-5 



+ 1.5E-2 
-7.6E-3 
+6.9E-3 
-9.8E-3 



+ 1.6E-1 
-1.9E-3 

+0.25 
-0.14 
+ 1.5E-3 
-4.7E-3 



1.2 
3.9 
1.1 
2.9 

0.11 
5.4E-4 

0.30 
3.4E-2 

0.30 
0.00 

2.34 

1.2E-2 
1348.1 
1231 



+0.003 
-0.02 

0.6 

-0.5 
0.04 
-0.03 

0.3 

-0.2 



+8 


3E- 


3 


-9 


.7E- 


4 


+6 


0E- 


5 


-7 


1E- 


5 


+2 


8E- 


2 


-1 


.5-2 




+2 


5E- 


3 


-1 


0E- 


2 



+3.9E-5 
-0.0 

+0.35 
-0.28 
+2.1E-3 
-1.1E-3 



0.7 

3.8 
0.7 
2.8 

0.089 
6.3E-4 

< 0.35 
2.3E - 3 

< 0.35 
5.8E-2 

1.61 

8.2E-3 
352.4 
300 



+4. 


7E- 


3 


-9. 


2E- 


3 


+3. 


6E- 


5 


-6. 


8E- 


5 



+3.7E-4 
-1.5E-3 



+ 1.9E-2 
-8.7E-3 



+ 1.6E-3 
-2.6E-3 



